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Abstract: 

We calculate the speed of sound Cg in an ideal gas of resonances, whose mass spectrum is 
assumed to have the Hagedorn form p{m) ~ m~"exp{6m}, which leads to singular behavior 
at the critical temperature = 1/b. With a = 4 the pressure and the energy density remain 
finite at Tc, while the specific heat diverges there. As a function of the temperature, the 
corresponding speed of sound initially increases similarly to that of an ideal pion gas, until near 
Tc resonance effects dominate, causing Cg to vanish as (Tc — T)^/^. In order to compare this 
result to the physical resonance gas models, we introduce an upper cut-off M in the resonance 
mass integration. Although the truncated form still decreases somewhat in the region around 
Tc, the actual critical behavior in these models is no longer present. 



1 Introduction 



The abundant formation of resonances of increasing mass and rotational degrees of freedom is 
one of the most striking features of strong interaction physics, which has attracted intense the- 
oretical attention for half a century or more. Even before the quark infrastructure of hadrons 
was known, a self-similar composition scheme, the statistical bootstrap model (SBM), led to an 
exponentially increasing resonance spectrum [1]. Shortly thereafter the dual resonance model 
(DRM) provided a description of hadron interactions in terms of the Regge resonance struc- 
ture in the different kinematic channels [2], which again produced an exponentially increasing 
spectrum [3]. In both these cases the basic feature is the underlying partition structure [4-6], 
leading to the form 



for the number of states (degeneracy) of resonances of mass m. This situation posed for the 
theory the two basic questions: what is the origin of so many resonances, and what is their 
effect on strong interaction thermodynamics? 

The first question was answered by the advent of the quark model and of quantum chromo- 
dynamics (QCD), whereby the resonances are the different possible excitation states of the 
(initially) three quark flavors, together with their antiquarks. The second one was answered 
by R. Hagedorn [1], who showed that an exponentially rising resonance mass spectrum leads 
to an upper bound on the temperature of hadronic matter, Tc = 1/6, with b as specified in 
the eq. ([1]). Beyond a certain point an increase in the energy density of the medium does not 
increase its temperature; instead, more and higher mass resonances are formed. The behavior 
of the medium at the "Hagedorn" temperature is determined by the power a in eq. ([1]). 
In particular, it is shown that for a > 4, the energy density itself remains finite at T^, only 
higher derivatives diverge there. The general structure of the resulting thermodynamics near 
Tc is thus that of critical behavior, associated with a phase transition leading to a new state of 
matter, which we today take to be the quark-gluon plasma (QGP) [7]. 

The aim of this work is to determine in some detail the onset of the critical behavior in strongly 
interacting matter, when coming from the confined hadronic side. The system here is an in- 
teracting hadron system; but it is known that if the interaction is resonance dominated, the 
medium of interacting elementary particles can be replaced by an ideal gas of all possible 
resonances [8,9]. If we consider systems of vanishing overall baryon number density, the as- 
sumption of resonance dominance appears to be well satisfied; only for high baryon density do 
non- resonant repulsive forces come into play [10]. 

We shall therefore consider here an ideal resonance gas of vanishing overall baryon density, and 
study in particular the temperature dependence of the speed of sound, Cs(T), in such a medium. 
By definition. 
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where s = {dP/dT)v is the entropy density and 
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is the specific heat at constant volume. If the resonance gas has a suitable degeneracy structure 
of the Hagedorn/DRM type, the specific heat diverges at the limiting temperature Tc, while 
the entropy density remains finite there. As a result, the speed of sound will vanish at Tc. It is 
thus expected that Cs(T) is at low temperature T that of an ideal gas of ground state hadrons 
("pions"), until T becomes high enough to bring resonance effects into play and eventually 
drive the system towards critical behavior. 

The speed of sound has long been considered as a sensitive indicator of critical behavior in 
strongly interacting matter. Early lattice studies indeed found the expected sharp dip of 
in the critical region [11,12]. Later, more precise studies seem to indicate a weaker dip [13], 
in accord with the expectation that the deconfinement transition is a rapid cross-over, rather 
than a genuine phase transition with singular behavior of thermal observables. The success of 
physical resonance gas models in the prediction of hadron abundances [14] provided the stimulus 
to also calculate the speed of sound in such schemes, using a grand canonical partition function 
of an ideal gas of all experimentally observed states up to a certain large mass. Again one finds 
a decrease of c^(T) as the temperature is brought near the expected critical point [13,15]. 

In the next section we will evaluate the thermodynamical functions, first for an ideal gas of 
pions and then for an ideal resonance gas with an exponentially growing mass spectrum of 
the form ([T]). We shall refer to such a system as a Hagedorn resonance gas. Following this, 
we calculate the resulting speed of sound and study in detail how is approaches the expected 
critical behavior. In the final section, we compare our results to those obtained for the physical 
resonance gas based on the actually observed resonances up to a mass limit of 2.5 GeV [16], as 
well as to the results from studies in finite temperature lattice QCD. 

2 Thermodynamics of a Hagedorn Gas 



For an ideal Boltzmann gas of identical scalar particles of mass mo and three charge states 
("pions") contained in a volume V, the grand partition function is defined as 



N 



This expression can be evaluated, giving 

\nZ{T,V) = 3 



VTmo'^ 
27r2 



^2(mo/T), 
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where K2{x) is the modified Bessel function of the second kind. From this equation we obtain 

'd\nZ\ Smo'^T^ 



for the pressure and 
BoiT) = 



Po{T) 



d\nZ 



dV 



T 



-K2{mo/T), 



dT )y 27r2 



for the energy density. Using these relations, 

£o(T)+Po(T) 3m2T 



27r2 
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4 K,{mo/T) + (^) K,{m,/T) 



(6) 



(7) 



(8) 



3 



provides the entropy density. 

We now extend this to an ideal Boltzmann gas of resonances, described by an exponentially 
increasing mass spectrum of the Hagedorn form ([1]). For the full spectrum of hadrons, including 
ground state pions as well the resonances, we thus write 

p(m) = 35(m — mo) + Am"'^ exp{m/Tc} 6'(m — 2mo). (9) 

Here the constant A (of dimension m^) provides the normalization of the resonance contribu- 
tions relative to that of the pions, and we have chosen the power term a = 4 in order to obtain 
second order singular behavior, as will be shown shortly. The ^-function assures that the reso- 
nance spectrum starts above the two-pion threshold. With the spectrum ([9]), the logarithm of 
the grand canonical partition function is now given by 

lnZ(T,V) = YI_ S^3mlK2{mo/T) + aJ dm exp{m/Tj ir2(m/T)| ; (10) 



the pressure thus becomes 



P(T) = Po(T) + / ^'^ exp{m/Tj K2{m/T), (11) 

with the pion gas pressure Po{T) given by eq. ([6]). Similarly, we obtain 

e{T) = SoiT) + dmm-^ exp{m/T,} 3 K2im/T) + ( ) Ki{m/T) , (12) 

for the energy density and 

AT r /m\ i 

s{T) = So{T) + — / dmm-' exp{m/rj A K2{m/T) + i-\Ki{m/T) . (13) 

for the entropy density; here So{T) and so(T) are the corresponding pion gas quantities and 

m- 

The Bessel functions appearing in these relations have the asymptotic form for large argument 



K^{m/T) = J^ exp{-m/T} [1 + 0(T/m)] . (14) 
V / m 

This shows that the integrals determining pressure, energy and entropy density all converge for 
for T < Tc, while diverging for T > T^. The convergence of the thermodynamics potentials at 
T = Tc is a direct consequence of our choice a = 4 in the resonance spectrum; for a < 7/2, 
both energy and entropy density would diverge there. In our case, one further temperature 
derivative is needed to cause a divergence, leading to a diverging specific heat at Tc. We shall 
return to the detailed form shortly. 

Given the analytical form ([9]) of the resonance spectrum, both eo{T) and e{T) can be evaluated 
numerically. For this, the resonance weight A has to be specified. With £{T^/T^ ~ 5, as 
obtained in lattice studies for the case of two quark flavors [17], we obtain A = 15.35 Tj^. Using 
this value, we get the energy density behavior of the ideal resonance gas shown in Fig. [1], where 
it is compared to that of the ideal pion gas. We see that around T/Tc ~ 0.6, resonances come 
significantly into play, so that e{T) begins to increase above the pion gas value, until at T^, 
resonances provide the dominant part of the energy density. 
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Figure 1: Energy density vs. temperature for an ideal pion gas (solid lower curve) and an ideal 
resonance gas (dashed upper curve) 

3 The Speed of Sound in a Hagedorn Gas 

The speed of sound Cg at constant volume is defined by eq. ([2]). For the case of the ideal pion 
gas we obtain from the eqs. ([3]) and f|T2|) 



C°(T) = 3so(T) + 
Together with the eqs. ([2]) and ^ this leads to 



K^inio/T). 
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mo^K2(mo/r) 



4T^K2{mo/T) + m,TK,{m,lT) 



(15) 



(16) 



for the ideal pion gas. Making use of the large argument limit f[T^ of the Bessel functions, we 
find for T ^ 0, the squared speed of sound vanishes linearly with T. The small argument 
limit 

Kn{x) = ^^• [l-O(x^)], (17) 

shows that for T oo, we get 1/3. The overall behavior below is shown in Fig. [2l 
For the Hagedorn gas, the same line of argument gives 

1 



Cv{T) = 3s (T) + 



27T^T 



/oo 
dm exp{m/Tc} K2{'m/T) 
2mo 



for the specific heat. As T — > Tc, the entropy density remains finite. However, the integral in 
eq. ( |T8i) diverges so that at the Hagedorn temperature we have the critical behavior given as 

Cv(T)~(T,-T)-°, (19) 

whereby the resulting critical exponent was found to be a = 1/2 [5]. 

An immediate further consequence of the diverging specific heat is that now the speed of sound 
vanishes at T^. To obtain the functional form of this behavior, we make use of eqs. (1131) . ([2]) 
and ([18]) to get 
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Figure 2: Speed of sound vs. temperature for an ideal pion gas (solid upper curve) and a 
Hagedorn gas (dashed lower curve) 

for the behavior in the Hagedorn gas. In this relation, the second term in the numerator 
diverges as T — >• T^, which in turn causes the speed of sound to vanish there. 

Making use of the value of A obtained above, we show in Fig. [2] the Hagedorn gas result 
compared to the ideal pion gas form. It is seen that the Hagedorn gas follows the pion gas 
form closely until about T/Tc ~ 0.2; from then on, any further increase of the average energy 
density is mainly turned into resonance masses, keeping the speed of sound roughly constant 
but non-zero. Only for temperatures very close to the critical point, for T/Tc ^ 0.9, do we reach 
a regime in which more and more of the energy goes into forming massive resonances. The 
truly critical region of the Hagedorn gas is thus restricted to a rather narrow band below the 
critical temperature. Magnifying this region, we show in Fig. [3] the behavior of as function 
of the expected critical form (T^ — T)^/^. It is seen that only in the last 1 - 2 % below the 
predicted behavior actually occurs. 
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Figure 3: Speed of sound of a Hagedorn gas (dashed lower curve) in the critical region, compared 
to the predicted critical behavior (Tc — T)^/^ (solid upper curve) 

We can obtain further information on the extent of the critical region by making use of the 
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interaction measure 



A(T) 



e(T) - 3P(T) 



(21) 



formed from the trace of the energy-momentum tensor; it vanishes for an ideal gas of massless 
constituents. In Fig. |U we show the modifications coming in for a resonance gas. Again we 
notice a rather large region of "normal" resonance effects, until very close to Tc, the critical rise 
sets in leading to a diverging slope of A(T) as T T^. 
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Figure 4: Interaction measure for an ideal pion gas (bottom curve) compared to that of a 
Hagedorn gas (top curve) 



4 The Low Mass Resonance Spectrum 

In this section we want to consider in some more detail the effect of the resonance spectrum 
form in the low mass range on the speed of sound Cs{T). Our choice of the power factor of 
p{m) in eq. ([1]) was dictated by the desired behavior of the thermodynamic observables at the 
divergence point T^. When a = 4, we had a finite energy density and a singular specific heat. 
The crucial feature for this was, however, the behavior at large m, for which a change in the 
threshold form at m = 2mo would not modify that property. The functional form for the 
spectrum ([1]) is shown in Fig. [5l It has a minimum at m = 4/6, before it begins its exponential 
rise. Thus it is evident that this minimum has no real physical meaning and can be replaced by 
a more reasonable threshold behavior of the resonance contributions. Therefore, we can replace 
p(m) in ([T]) by 

2mo' 



Ps 



m] 



1 - 



m 



m 



exp{6m} 



(22) 



with some integer constant power s > 0. Clearly Ps{m) still retains the same critical behavior. 
The spectral form produced by ( !22l) for s = 2 is included in Fig. [5l Since the region of the 
low mass resonances contributes very little to the functional form of the speed of sound, the 
changes that such a threshold factor produces in cl(T) are effectively negligible. 
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Figure 5: Resonance spectrum for eq. ([T]) with a = 4, compared to the form given by eq. (|22ll . 

5 The Physical Resonance Gas 

We now want to compare our results to those obtained in studies of an ideal (Boltzmann) gas 
of actual physical resonances. The relevant partition function thus becomes the sum over all 
states listed by the Particle Data Compilation [16] up to some mass value of M, above which 
the experimental information is too incomplete to be of use. Typically one has M = 2 — 3 GeV. 
For each state of mass rrii, we include a degeneracy gi determined by the spin, isospin, baryon 
number and strangeness degrees of freedom. The partition function is then given by 

In Zp{T,V) = Y,ln Z},{T,V), (23) 

i 

with 

In ZUT, ^) = ^ / d'pexp |-^^^^^| = ^ m^iT^K/T), (24) 
specifying the contribution of a resonance of rrii. 

Before we calculate the speed of sound for such a "physical" resonance gas, let us see what effect 
a mass truncation in the resonance integration has for the Hagedorn gas. We thus consider the 
form obtained if we replace in eq. (fTOj) the integral over [2mo, oo] by one over [2mo,M]. The 
results are shown in Fig.[6]for several values of M. Since the truncated system has no divergence 
at T = Tc, that is since the specific heat now remains finite there, we can continue the curves 
also to values above T^. Eventually, for any finite M, they will converge to = 1/3. 

We now turn to the physical resonance gas and calculate the speed of sound for several values 
of the upper limit M in resonance mass. The results are shown in Fig. [71 Qualitatively, they 
agree quite well to what was found for the truncated Hagedorn gas. It is seen that in both 
scenarios, the inclusion of further heavy resonances somewhat lowers the speed of sound in the 
temperature range < T < 200 MeV, but it never produces a sharp minimum. 

Futher inspection does show, however, a difference in the behavior of the two formulations. In 
Fig. El we show the speed of sound at T = Tc in the Hagedorn gas for increasing resonance 
mass bounds. It is seen to converge to zero for M — oo, which means that it is only the mass 
bound which prevents critical behavior. 

In the case of the physical resonance gas, we also have to consider the form of the degeneracy. 
In the degeneracy factor Qi, only the spin variable Jj can result in an unbounded increase with 
i. For resonances on linear Regge trajectories, we have 

Ji ~ a'mf, (25) 
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Figure 6: Speed of sound in a Hagedorn gas with resonance mass truncation, m < M. 
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Figure 7: Speed of sound in a physical resonance gas with upper hmits M in resonance mass. 



where a' ~ 1 GeV is the Regge resonance slope. Hence the spin degeneracy (2 Jj + 1) implies 



g{mi) ~ m. 



(26) 



It is clear that with such a degeneracy, the partition function will remain analytic for all values 
of T, so that we will never encounter critical behavior, no matter what the upper bound on the 
mass is. This explains why, as seen in Fig[9l an increase of M appears to produce a convergent 
finite value of cj.(Tc). 

Thus it is not only the upper limit in resonance mass which prevents the physical resonance gas 
from ever showing critical behavior. The crucial factor is the use of a degeneracy determined 
by quantum number degrees of freedom only. Such a degeneracy grows, as we saw above, as 
m^, and as a consequence, the partition function always remains analytic. The critical behavior 
arising for a Hagedorn resonance gas is due to a degeneracy growing exponentially in m. This, 
in turn, is a consequence of counting as distinct degenerate states all partitions of a given 
resonance into all other possible resonances. As we noted above, both the SBM and the DRM 
obtain the exponential mass degeneracy from the corresponding partition problem. 
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Figure 8: Speed of sound for a Hagedorn gas as function of the upper limit M in resonance 
mass. 
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Figure 9: Speed of sound at T = 175 MeV in a physical resonance gas, as function of upper 
resonance mass limit M. 

Do these considerations have any physical implications? They clearly would, if deconfinement 
were an actual phase transition with a diverging specific heat. If it is only a rapid cross-over, 
or if the transition should lead to a finite specific heat at Tc (such as produced in the 0(4) 
universality class), the behavior of the speed of sound obtained from the physical resonance gas 
even at can remain a good approximation. 

On the other hand, we have also seen that the actual critical region in the Hagedorn gas is a 
very narrow band just below Tc. Up to a temperature of T ~ 0.95 Tc, the observed behavior is 
due to the normal resonance pattern, and here physical resonance gas and Hagedorn gas can 
certainly agree. This leaves us with an amusing question: is the temperature obtained in abun- 
dance analyses of hadron species indeed the temperature of the confinement /deconfinement 
transition, or does the replacement of an interacting hadron gas by a physical resonance gas 
break down very close to a critical point? 
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